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Summary. - The connection between the strictly isospectral construction in 
supersymmetric quantum mechanics and the general zero mode solutions of the 
Schroedinger equation is explained by introducing slightly generalized first-order 
intertwining operators. We also present a multiple-parameter generalization of 
the strictly isospectral construction in the same perspective. 
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1. - Consider the second order, one-dimensional, differential operator of 
Schroedinger type 

H- = -D 2 + V-(x), (1) 

where D = 4- and V-(x) is a "bosonic" type potential. We are interested in 
solving the eigenvalue equation associated with H_ by means of factoring it as 
the product of two first order operators. If we consider the ground state case, 
and moreover take the ground state energy as zero, we get the second order 
linear differential equation = 0. The general solution of such a linear 

equation will be a linear composition of the two linearly independent particular 
solutions, usually known as zero modes in the literature. We shall suppose 
that at least one of the particular solutions, say u, is nodeless and normalizable 
(a true zero mode), a case denoted as unbroken supersymmetry If u is 
supposed to be known, it allows us to factorize H- in the following way Q. 
Defining W = —Inu as the bosonic superpotential one gets 

H- = -D 2 + V-(x) = A^A = (-D + W')(D + W), (2) 

where W = ^ and V-{x) = —W" + W' 2 = The true zero mode u is 
annihilated by the A operator Au — 0. 

We can construct a second linearly independent particular solution t;_ from 
the first one u, namely B 

r x l 

v. = uj - 2 . (3) 

This solution is not annihilated by the A operator, Av- ^ but acting on 
Av- — — is zero, since v_ is a zero mode of H_. 

The general solution annihilated by ii_ can be constructed as a linear com- 
bination of u and i>_ as $_ = au + bv- = b(^u + u_) cx (X s u + The range 
of variation of A s can be extended such that A s G (-co, oo); A s = selects the 
u_ solution and X s ±oo the u one. If we factorize H— by the usual proce- 
dure but now using the general solution $_(A S ) to construct a A s -dependent 
superpotential W(X S ) — — ln<f>_(A s ), the factoring operators A and A^ will also 
depend on A s , but nevertheless the A s dependence disappears in the final result 
and one ends up with the original potential V- 

H_(\ s ) = (-£> + W'(X S )){D + W'(X S )) = A\X S )A(X S ) = -D 2 + V-(x). (4) 

The supersymmetric "fermionic" partner of (|l|) is constructed by interchang- 
ing the order of the factors 0] 

H+ = AA f = {D + W')(-D + W') = -D 2 + W' 2 + W" = -D 2 + V+(x). (5) 

The "fermionic" potential differs from the "bosonic" one by — 2W", and, as is 
well known, both operators H + and iJ_ have the same spectrum apart from the 
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ground state which disappears in as a matter of fact, one can write down 



H_\n >= E n \n >^ A^A\n >= E n \n A(A*A\n >) = A(E n \n >) 
=> AA\A\n >) = E n (A\n >) => H+\n >= E n \n >, (6) 

and for the ground state, ^4|0 >= 0. 

One particular solution of the differential equation H + Q + — is ^, which, 
as said before, is annihilated by A^; we can construct the second linearly inde- 
pendent solution by the same procedure as for H- 

1 r x 

,+ = - u\ (7) 



u „ 

and the general solution of H + ^> + = reads 

<MA S ) = A s i + v + = ~ f\ s + J X u 2 ^j , (8) 

which we call a one-parameter fermionic zero mode. 

Now again, if we use the general solution (|^) to factorize the fermionic Hamil- 
tonian we can do it in an infinite number of ways, depending on the value of 
X s but again as in the bosonic case all factorizations are producing only the 
original potential. 

The interesting property of <&+ is that its inverse has the form of the gen- 
eral wave function which is used to generate one-parameter families of strictly 
isospectral potentials to the original bosonic one in the supersymmetric ap- 
proach ||. Thus 

can be used to construct a family of Hamiltonians 

H(X S ) = A\X S )A(X S ) = -D 2 + V-(\ B ), (10) 

where A(X S ) = D + W(X S ) and W(X S ) = -ln(^ A J- These Hamiltonians are 
strictly isospectral (there is absolutely no difference in their energy spectra, 
though there are differences in their wave functions) since they can be con- 
structed by means of the general solution of the Riccati equation; the isospec- 
tral family includes the original potential V_ corresponding to the superposition 
quotient X s — > ±oo. We briefly present this method in the following. 

2. - Consider the "fermionic" Riccati equation (FRE) y = — y 2 + f(x) 
with the known particular solution yo and let y\ = wi + yo be the general 
solution. By substituting y\ in FRE one gets the Bernoulli equation —w 1 = 
w 2 + (2yo)w\. Furthermore, using W2 = 1/wi, the simple first-order linear 
differential equation w 2 — (2y n )w 2 — 1 = is obtained, which can be solved by 
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employing the integration factor / = e / 2yo , leading to the solution wi = 

+ X r ), where X r occurs as an integration constant. Thus, 
the general solution of FRE is 

e-f 2yo f 

2/1 = 3/0 H ft— =2/o + - — -7x777 ■ (11) 

X r + J x e-J 2yo X r + I (/) 

One can see easily that the particular Riccati solution tjq corresponds to Witten's 
superpotential jlj] while the general Riccati solution y\ — yo + - -^-ps - is of 

Mielnik type ||. Also u = f 1 ^ 2 can be normalized being the ground state 
wavefunction (the true zero mode), and — 2y (= — In/ 1 / 2 ) is the Darboux 
transform part of the Schroedinger potential. Moreover, one can see easily that 
the modes 

*a„ = J r = - r 2 (12) 
A,. + J / A r + J u 2 

can be normalized and therefore are the ground state wavefunctions (true zero 
modes of the bosonic family) corresponding to Mielnik's parametric superpo- 
tential. Also, —2y 1 can be thought of as the general Darboux transform part in 
the potential generating the bosonic strictly isospectral family, which reads 

This family of potentials can be seen as a continuous deformation of the original 
potential, because this is included in the infinite limit of the deforming param- 
eter. All these relationships are supplementary material to the core of the 
"entanglement" between Riccati and Schroedinger equations, which has been 
recently emphasized by Haley Q. 

If now, the two A constants we have used are identified, i.e., A r = A s = A, 
then of course ^>\ r = . Since X s is a superposition quotient while A r is an 
integration constant, and both can be fixed through boundary conditions, we 
can see that the identification is sound. This also gives a clue on the connection 
between Schroedinger and Riccati methods at the level of their general solutions. 

We have now to answer the important question why the factorization based 
on 4" \ r gives a nontrivial and therefore a physically relevant case and what 
mathematical feature lies behind (^). These issues are addressed in the follow- 
ing. 

3. - The mathematical background which is needed here belongs to the in- 
tertwining operator transformations studied mathematically by Moutard, Dar- 
boux, Ince, Delsart, Lions, and others. For a detailed analysis see ||. Inter- 
twining has been used by Pursey in his studies of the combined procedures for 
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generating families of strictly isospectral Hamiltonians Q , whereas Anderson JtJ 
applied matrix intertwining relations to the Dirac equation showing that their 
structure is described by an N — 4 superalgebra. 

Two operators Lq and L\ are said to be intertwined by an operator T if 

L X T = TL . (14) 

If the eigenfunctions ip a of L are known, then from the intertwining relation 
one can show that the (unnormalized) eigenfunctions of L\ are given by ip% = 
TifQ. The main problem in the intertwining transformations is to construct 
the transformation operator T. One-dimensional quantum mechanics is one 
of the simplest examples of intertwining relations since Witten's transformation 
operator T qm = T\ is just a first spatial derivative plus a differentiable coordinate 
function (the superpotential) that should be a logarithmic derivative of the 
true bosonic zero mode (if it exists), but of course higher-order transformation 
operators can be constructed without much difficulty j^] . 

Thus, within the realm of one-dimensional quantum mechanics, writing 7\ = 

D — — , where u is a true bosonic zero mode, one can infer that the adjoint 

operator t\ = —D — 2- intertwines in the opposite direction, taking solutions 
of L\ to those of Lq 

ipo = T x Vi ■ (15) 

In particular, within quantum mechanics, Lq = H- and L\ = H + and although 
the true zero mode of H_ is annihilated by T\ , the corresponding (unnormalized) 
eigenfunction of H + can nevertheless be obtained by applying T\ to the other 
independent zero energy solution of H— . It is now not very difficult to see that 
a slightly modified transformation operator of the adjoint type 

T- = -D + *l - |A (16) 

when applied on the general zero mode solution 4> + of H + will produce precisely 
the strictly isospectral family of bosonic zero modes (T A ~ ( E ) + = ^\). On the 
other hand, the corresponding transformation operator sending general bosonic 
zero modes to general fermionic ones reads 



7T ^ (17) 



Moreover, one can start with one of the strictly isospectral bosonic zero 
modes = x , = u i (i- e -> fixing A = Ai) and repeat the strictly isospec- 
tral construction, or in other words, using a new two-parameter (of which only 
the second one is a free parameter) fermionic zero mode of a Hamiltonian H ++ 

$+(Ai,A 2 ) = ^-(a 2 + J u\\ , (18) 
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whose inverse ^ r Ai,A 2 = x+J* i = M2 ' ^ 2 e l" 00 ; 00 ); ^ s J us ^ a two-parameter 
family of bosonic zero modes. The resulting two-parameter family of strictly 
isospectral potentials will be 

Vm M = V Xl (x) - 2 2 - + 2 , (19) 

A 2 +J ui, (A 2 + J uj) 2 

where V\ 1 is given by (|l3|), with \ r = Ai. Generalized formulas for multiple- 
parameter cases can be easily written down, as well as a multiple parameter 
form of the modified transform operators. 

As we mentioned, the general bosonic zero modes, u\, m 2 , Us,...Ui, can 
be normalized and therefore turned into true zero modes. Their normalization 
constants are of the type Ni = y/Xi(Xi + 1), where i = 1,2,3,..., whenever 
one is performing the normalization at the previous (i — l)-parameter step and 
considering u = u as a true zero mode (i.e., N — 1, and thus already normal- 
ized to unity). At an arbitrary z-level, the parameter dependence of Ui, when 
expressed in terms of only uq is a rather complicated type of denominator pos- 
sessing nesting integrals, which physically may be considered as a modulational 
factor endowed to the multiple-parameter, bosonic, true zero mode through 
the mathematical procedure. The parametric normalization deletes the interval 
[—1,0] from the parameter space; at the —1 limit, one can make a connec- 
tion with the Abraham- Moses isospectral technique ||, whereas at the limit 
the connection can be done with another isospectral construction developed by 
Pursey 0. Moreover, since the strictly isospectral supersymmetry obviously 
may introduce singularities in both wavefunction and potential, usually the ac- 
tive authors in the field are discarding those values of the deformation parameter 
for which those singularities are occuring. If |A s j ng | > 1 the excluded interval 
in the parameter space is [— A s i ng , A S i ng ], and therefore the connections with the 
Abraham-Moses and Pursey methods are lost in this case. Usually, at least 
one of the limits is lost, as for example in the harmonic oscillator case where 
the excluded interval is [— -^] implying the loss of the Pursey limit. Also, 
the strictly isospectral supersymmetric method (as well as the other isospectral 
procedures) may be considered as a way of allowing for some physical effects of 
the irregular Schroedinger (vacuum) solutions in quantum mechanics, which is 



not so much of a nonsense nowadays 1 10 . Indeed, the wrong vacuum is wrong 



only because of its asymptotic behaviour; in some cases, the strictly isospectral 
constructions show that the asymptotic behaviour of the wrong vacuum can be 
tamed. 

In conclusion, we presented in some detail the parameter dependence of the 
zero energy sector of the unbroken supersymmetric quantum mechanics, which 
in applications is usually fixed through some sort of boundary conditions. 
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